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ABSTRACT 



Aims. We have been carrying out a precise radial velocity (RV) survey for K giants to search for and study the origin of the low- 
amplitude and long-periodic RV variations. 

Methods. We present high-resolution RV measurements of the K2 giant HD 66141 from December 2003 to January 2011 using the 
fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at Bohyunsan Optical Astronomy Observatory (BOAO). 
Results. We find that the RV measurements for HD 66141 exhibit a periodic variation of 480.5 ± 0.5 days with a semi-amplitude 
of 146.2 ± 2.7 m s -1 . The Hipparcos photometry and bisector velocity span (BVS) do not show any obvious correlations with RV 
variations. We find indeed 706.4 ± 35.0 day variations in equivalent width (EW) measurements of H ff line and 703.0 ± 39.4 day 
variations in a space-born measurements 1.25yu flux of HD 66141 measured during COBE/DIRBE experiment. We reveal that a mean 
value of long-period variations is about 705 ±53 days and the origin is a rotation period of the star and variability that is caused 
by surface inhomogeneities. For the 480 day periods of RV variations an orbital motion is the most likely explanation. Assuming a 
stellar mass of 1.1 ±0.1 M e for HD 66141, we obtain a minimum mass for the planetary companion of 6.0 ± 0.3 M Jup with an orbital 
semi-major axis of 1.2 ± 0.1 AU and an eccentricity of 0.07 ± 0.03. 

Key words, stars: planetary systems - stars: individual: HD 66141 - stars: giant - technique: radial velocity 



1. Introduction 

So far, more than 570 exoplanets have been detected and roughly 
80% of them have been found via high resolution spectroscopic 
measurement of radial velocity (RV) of spectral lines. Besides 
RV measurements, high resolution spectroscopy simultaneously 
offers various information on the chemical composition of host- 
ing stars. Among the planetary companions discovered by the 
RV method, only ~ 30 exoplanets have been detected around 
giant stars. Cool evolved stars such as K giants are suitable can- 
didates for precise RV measurements because they have sharp 
and sufficient spectral lines compared to those of early spec- 
tral class intermediate-mass and rapidly-rotating main-sequence 
(MS) stars. As the star moves towards giants with convective 
outer envelopes, the role of surface inhomogeneities (chromo- 
spheric activities or stellar spot modulations) increase. Thus, it 
is difficult to distinguish the RV variation due to a planetary com- 
panion from that due to an inhomogeneity in evolved stars and 
therefore, giants still remain an undeveloped region in exoplanet 
surveys. The suggestion of planetary companions in K giants 
(Hatzes & Cochran 1993) was first confirmed by Hatzes et al. 
(2005, 2006) and Refferet et al. (2006) after 13 years of obser- 
vation. 

Currently, a number of groups are conducting exoplanet sur- 
veys around giants (Frink et al. 2002; Sato et al. 2003; Setiawan 
et al. 2003; Hatzes et al. 2005; Johnson et al. 2007; Lovis & 
Mayor 2007; Niedzielski et al. 2007; Dollinger et al. 2009; Han 
et al. 2010). For the past eight years, we have conducted precise 
RV measurements of 55 K giants. Here, we present a long-period 



and low-amplitude RV variation of a K giant HD 66141. Data 
observations and analysis are presented in Sect. 2. In Sect. 3 and 
Sect. 4, we describe the properties of HD 66141 and analyze the 
period search, respectively. The nature of the RV variations is 
investigated in Sect. 5. Finally, we discuss the work presented in 
this paper in Sect. 6. 

2. Observations and analysis 

We started a precise RV survey using 1.8-m telescope at BOAO 
in 2003 to search for exoplanets and to study the asteroseismol- 
ogy of 55 K giants. We have reported two new exoplanets (Han 
et al. 2010; Lee et al. 201 1), a confirmation of an exoplanet (Han 
et al. 2008), and an oscillating star (Kim et al. 2006). 

We acquired 54 spectra for HD 66141 (= HR 3145 = 
HIP 39311) from December 2003 to January 2011 using the 
fiber-fed high-resolution (R = 90 000) Bohyunsan Observatory 
Echelle Spectrograph (BOES; Kim et al. 2007) attached to the 
1.8-m telescope at Bohyunsan Optical Astronomy Observatory 
(BOAO) in Korea. An iodine absorption cell (L) was used to 
provide the precise RV measurements. Each estimated signal-to- 
noise ratio (S/N) at the L wavelength region is about 200 - 250, 
with typical exposure time ranging between 240 and 480 sec- 
onds. The RV measurements forHD 66141 are listed in Table|2] 

The reduction was carried out using IRAF (Tody 1986) soft- 
ware and DECH codes (Galazutdinov 1992). IRAF was used for 
image processing and the extraction of 1-D spectra, and a con- 
tinuum process was performed using the DECH code. L anal- 
yses and precise RV measurements were undertaken using the 
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Table 1. Stellar parameters for HD 66141. 



Parameter 


Value 


Reference 


Spectral type 


K2III 


Perryman et al.(1997) 


m v (mag) 


4.39 


Perryman et al.(1997) 


M v (mag) 


-0.15 


Perryman et al.(1997) 


B-V(mag) 


1.26 


Massarotti (2008) 


age (Gyr) 


6.84 ± 1.39" 


Derived 


Distance (pc) 


80.9 ± 6.3 


Famaey et al. (2005) 




80 


Massarotti et al. (2008) 


RV(kms-') 


71.57 ± 0.01 


Famaey et al. (2005) 


Parallax (mas) 


12.49 ± 0.98 


Perryman et al.(1997) 




12.48 


Kharchenko et al. (2007) 




12.84 ± 0.25 


van Leewen (2007) 


Diameter (mas) 


2.81 ± 0.15 


Richichi & Percheron (2002) 


Teff (K) 


4305 


Massarotti et al. (2008) 




4323 ± 15 


This work 


[Fe/H] 


-0.36 


Massarotti et al. (2008) 




-0.323 ± 0.034 


This work 




2.3 


Massarotti et al. (2008) 




1.92 ±0.07 


This work 




1.78 ± 0.04° 


Derived 




22 


Massarotti et al. (2008) 




21.4±0.6 fl 


Derived 


Af* (M G ) 


1.1 +0.1" 


Derived 


U (L ) 


174 


Massarotti et al. (2008) 


v rot sin i (km s -1 ) 


2.5 


Fekel (1997) 




1.1 ± 1.0 


de Medeiros & Mayor (1999) 




4.7 


Massarotti et al. (2008) 


P mt 1 sin i (days) 


253 - 938 


This work 


Vmicro (km S~') 


1.45 ±0.08 


This work 



Derived using an online tool (http://stevoapd.inaf.it/cgi-bin/param I. 



RVI2CELL (Han et al. 2007) developed at the Korea Astronomy 
& Space Science Institute (KASI). The RVI2CELL adopts basi- 
cally the same algorithm and procedures described by Butler et 
al. (1996). However, for modeling the instrument profile we used 
the matrix formula described by Endl et al. (2000). We solved 
the matrix equation using singular value decomposition instead 
of the maximum entropy method adopted by Endl et al. (2000). 
The code have been used in several cases (Kim et al. 2006; Han 
et al. 2008, 2010; Lee et al. 2008, 2011) 

The RV standard star r Ceti shows an rms scatter of 6.7 m 
s _1 over the time span of our observations and demonstrates the 
long-term stability of the BOES (Lee et al. 201 1). 

3. The properties of HD 66141 

HD 66141 is an IAU bright (V = 4.39, K2 III) RV standard star 
adopted by IAU Commission 30 (Pearce 1955). It has been ob- 
served very frequently since and multiple RVs of the object are 
available in the literature (e.g. RV = 71.6 + 0.3 m s _1 [Udry et 
al. 1999a]; 71.40 + 0.15 m s" 1 [Eaton & Williamson 2007]). The 
RV stability of HD 66141 and other IAU standards were indeed 
under doubts during last decades (Batten 1983). Over 20 years 
of its and other IAU standard observation with CORAVEL were 
summarized in Udry et al. (1999a). While, authors did not find a 
regular variability of RVs in HD 66141, Udry et al. (1999b) did 
not include it to a list of recommended standards for future use 
due to the fact that as a giant it exhibits some RV variability at 
higher precision level. 

We determined the atmospheric parameters of HD 66141 
directly from our spectra. By using 266 measured equivalent 
widths (EWs) of Fe I and Fe II lines, we determined T e e, [Fe/H], 



Table 2. RV measurements for HD 66141 between December 
2003 and January 2011. 



JD 


ARV 


±<T 




JD 


ARV 


±<T 


-2 450 000 


m s 


m s 


-2 450 000 


ms~' 


ms~' 


2976.251556 


-49.8 


9.7 


3459. 


.015331 


-59.9 


7.3 


2977.315941 


-123.3 


8.1 


3459. 


.080221 


-74.1 


7.7 


2977.324622 


-123.6 


8.1 


3729. 


.107766 


189.1 


9.2 


2978.186977 


-35.9 


9.1 


3759. 


.164025 


132.2 


7.2 


2978.205694 


-36.5 


9.0 


3778. 


.200411 


142.8 


8.2 


2980.409693 


-75.4 


14.6 


3779. 


.170158 


167.8 


7.4 


3045.039195 


-15.0 


6.7 


3818. 


.009412 


40.4 


6.7 


3045.047343 


-16.6 


6.7 


3819. 


.971947 


64.5 


7.7 


3046.104835 


-51.8 


9.0 


4036. 


.337614 


31.3 


6.5 


3046.115923 


-48.9 


7.9 


4123. 


.092626 


162.7 


7.2 


3046.123145 


-45.7 


11.2 


4126. 


.108709 


172.3 


6.5 


3048.085051 


0.8 


7.3 


4214. 


.007594 


203.0 


10.8 


3048.092203 


8.9 


8.3 


4396. 


.340192 


-125.4 


6.7 


3072.070064 


5.0 


10.2 


4452. 


.360964 


-146.5 


7.3 


3072.078211 


16.6 


10.9 


4506. 


.168231 


-79.0 


7.2 


3095.998435 


16.9 


7.1 


4536. 


.065764 


-4.7 


6.7 


3132.993871 


76.4 


7.2 


4538. 


.042117 


18.3 


6.5 


3332.318368 


-33.1 


7.3 


4755. 


.368868 


27.4 


7.4 


3332.325255 


-35.6 


6.9 


4847. 


220274 


-134.2 


6.9 


3354.234041 


34.0 


6.9 


4880. 


.152070 


-144.2 


8.5 


3395.246936 


-53.8 


11.9 


4929. 


020161 


-142.7 


7.1 


3430.046636 


-21.6 


7.3 


5171 


.229126 


174.3 


8.2 


3432.981185 


-40.8 


9.2 


5248. 


.098927 


27.2 


6.4 


3432.991011 


-44.8 


8.1 


5251. 


.084871 


16.6 


6.6 


3433.115135 


-63.5 


8.0 


5456. 


.347802 


-4.0 


8.4 


3433.126349 


-58.0 


7.6 


5554. 


.319097 


181.7 


15.1 


3433.183059 


-68.1 


8.8 


5581. 


.132641 


144.7 


8.6 



log g, and v m i clo of the star using the program TGVIT (Takeda et 
al. 2005). To estimate the stellar radius, mass and surface gravity, 
we used an online tool (http://stev.oapd.inaf.it/cgi-bin/param), 
which is based on theoretical isochrones (Girardi et al. 2000; 
J0rgensen & Lindegren 2005; da Silva et al. 2006). The result is 
R+ = 21.4 + 0.6 R , log g = 1.78 + 0.04, and M* = 1.1 + 0.1 

The v rot sin i values of HD 66141 are from de Medeiros 
& Mayor(1999), Massarotti et al. (2008) and Fekel (1977). 
Two former author measured the rotational velocities using the 
spectroscopic line broadening methods and the cross-correlation 
spectrometers. Fekel (1977) used a high-resolution spectra. 
These determinations were also corrected for macroturbulence 
and instrumental broadening. The rotational velocity of 2.5 km 
s _1 (Fekel 1997) is close to the 1.1 + 1.0 km s -1 (Medeiros 
& Mayor 1999) but different with determination of 4.7 km s _1 
(Massarotti et al. 2008). With these determinations we have a 
range of uncertainty of v ro t sin i measurements 1.1 — 4.7 km s _1 . 
Adopting a stellar radius of 21.44 + 0.58 R e and ignoring the sin 
i we get the uncertainty range for the upper limit of the rotation 
period 

P rot = 2;rfl* / (v rot sin = 253 - 938 days. 

The basic stellar parameters for HD 66141 are summarized in 
Tabled] 

4. Period search 

RV measurements for HD 66141 are presented in Figure Q] and 
have a standard deviation of 94.5 m s _I , which is 14 times larger 
than the RV standard star t Ceti. The long-term periodic vari- 
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Fig. 1. RV curve (top panel) and rms scatter of the residual (bot- 
tom panel) for HD 66141 from December 2003 to January 201 1. 
The solid line is the orbital solution with a period of 480.5 + 0.5 
days and an eccentricity of 0.07 + 0.03. 
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Fig. 2. The Lomb-Scargle periodogram of the RV measurements 
for HD 66141. The periodogram shows a significant power at 
a frequency of 0.002081 + 0.000002 c d -1 corresponding to a 
period of 480.5 + 0.5 days (top panel) and after subtracting the 
main frequency variations (bottom panel). The horizontal dotted 
lines indicate an FAP threshold of 1 xl0~ 3 (0.1%). 



ability in the data is indeed visible by eye. The Lomb-Scargle 
periodogram for unequally spaced data (Lomb 1976; Scargle 
1982) was applied to the RV time series of HD 66141 in or- 
der to find an accurate period. The frequency uncertainties are 
standard least-square uncertainties of a sine-wave fit to the data. 
The periodogram shown in Figure|2](top panel) denotes a signif- 
icant power at /j RV = 0.002081 + 0.000002 cd" 1 (Pf = 480.5 
+ 0.5 days). We determined the significance of the period by 
calculating the false alarm probability (FAP) for the dominant 
period by a bootstrap randomization technique (Kiirster et al. 
1999). We computed that among 200 000 trials the highest peak 
of periodograms, thus we found an FAP of less than 10 -6 for 
/. RV . Our criterion for statistically significant signal in the data is 
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Fig.3. The 1.25/v flux intensity measurements of HD 66141 
phased to the period of 256 (top panel) and 703 days (bottom 
panel). Epochs of maximum flux for 256-day and 703-day peri- 
ods are JD 2447945.3 and JD 2448098.0, respectively. 



that it must have FAP< 10~ 3 . We note, that the 480.5 day period 
is within the uncertainty range 253 - 938 days for the rotation 
period found in the previous Section. 

Figure [2] (bottom panel) shows the Lomb-Scargle peri- 
odogram of the residuals after removing the best fit; it exhibits 
no statistically significant peaks at the domain of interest. 



5. Origin of the RV variations 

Evolved stars exhibit pulsations as well as surface activity, and 
they result in low-amplitude RV variabilities on different time 
scales. While short-term (hours to days) RV variations have been 
known to be the result of stellar pulsations (Hatzes & Cochran 
1998), long-term (hundreds of days) RV variations with a low- 
amplitude may be caused by stellar pulsations, rotational mod- 
ulations by inhomogeneous surface features, or planetary com- 
panions. To establish the origin of the period for HD 66141, we 
examined 1 ) the Hipparcos and COBE/DIRBE infrared photom- 
etry, 2) the stellar chromospheric activity, 3) the spectral line 
bisectors, and 4) the orbital fit. 

5.1. Hipparcos and COBE/DIRBE infrared photometry 

We analyzed the Hipparcos photometry (ESA 1997) for HD 
66141 to search for any possible brightness variations due to the 
rotational modulation of stellar spots. For three years, between 
JD 2447966 and JD 2448754, the Hipparcos satellite obtained 
45 photometric measurements for HD 66141 and the data main- 
tained a photometric stability down to the rms scatter of 0.0078 
magnitude, corresponding to 0.17% variations. Figure [9] shows 
the Lomb-Scargle periodogram of these measurements. There is 
no significant peak near the frequency of / RV = 0.00208 c d _1 . 
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3955 3960 3965 3970 3975 3980 3985 
Wavelength (A) 

Fig. 4. Ca II H spectral region for HD 66141. It shows that the 
Ca II H core feature exhibits no emission at the line center. 
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Fig.5. Line profile near H a region for HD 66141. Two vertical 
dashed lines denote the range of H a to measure the EWs. 



Two statistically insignificant peaks are visible at frequencies 
/™ p = 0.00233 + 0.00002 c d" 1 and f* w = 0.00274 + 0.00002 c 
d -1 with FAPs of ~ 5 xl0~ 3 that most likely belongs to 1/365.25 
day = 0.00274 cd~' artifacts. So, we conclude that Hipparcos 
photometry does not resolve any signals of surface features like 
spot. 

Nearly at the same time, between JD 2447973 and 
JD 2449128, HD 66141 was measured in the near-infrared (NIR) 
1.25, 2.2, 3.5, and 4.9yU bands by NASA's COBE (Cosmic 
Backgroud Explorer) satellite using DIRBE (Diffuse Infrared 
Background Experiment) instrument. The total N = 68 weekly 
averaged fluxes in each band were extracted for HD 66141 from 
COBE/DIRBE archives (Price et al. 2010) which we used for 
analysis. The Lomb-Scargle periodogram analysis of 2.2, 3.5, 
and 4.9p fluxes does not reveal any significant signals in the 
domain of interest, while the 1.25yU shows the statistically sig- 
nificant signal (FAP < 10~ 5 ) at the frequency of /j IR = 0.00390 
± 0.00008 cd l (Pf = 256.1 ± 5.4 days) with semi-amplitude 
Al = + 4.0 Jy. The second statistically significant peak (FAP < 
10~ 4 ) is at the lower frequency of / 2 IR = 0.00142 ± 0.00008 c d _1 
(Pf = 703.0 ±39.4 days). 

After removal from the original data of the 0.00390 or 
0.00142 cd _1 signals the residuals do not show any significant 
variations (Figure|9]). Figure[3]shows the phase curve of the orig- 
inal data folded with two periods P IR = 256 days (upper panel) 
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Fig. 6. The examinations of the origin of the RV variations for 
HD 66141. JD vs. U a EW variations (top panel) and the BVS 
variations (bottom panel), respectively, from December 2003 to 
January 201 1. 



and P' = 703 days (bottom panel). We note that 256-day and 
703-day periods are within uncertainty limits for the rotation pe- 
riod of 253 - 938 days found in Section 3. 

5.2. Chromospheric activity 

The EW variations of Ca II H and H^ lines are frequently used as 
chromospheric activity indicators. The emissions in the Ca II H 
core are formed in the chromosphere and show a typical central 
reversal in the existence of chromospheric activity (Pasquini et 
al. 1988; Saar & Donahue 1997). The existence of extra emis- 
sion at the center of the line implies that the source function in 
the chromosphere is larger than in the photosphere. This rever- 
sal phenomenon is common in cool stars and is intimately con- 
nected to the existence of a convective envelope and magnetic 
activity. Unfortunately, the Ca II H line region for HD 66141 
does not have enough S/N to estimate EW variations (Figure SJ. 
However, it is large enough to check the emission feature in the 
Ca II H line core, and we did not find any emission. 

Because the H„ absorption originates in the upper layers of 
stellar atmosphere and is also sensitive to stellar activity (Kiirster 
et al. 2003), we used the H EW to measure the variations. 
Owing to sparse blending lines, a weak telluric line, and a nar- 
row H,, absorption line (Figure |5j, it is easy to estimate the H a 
EW. We measured the EW using a band pass of + 1 .0 A centered 
on the core of the H„ line to avoid nearby blending lines (i.e. Ti I 
6561.3, Nail 6563.9, and ATM H z O 6564.0 A). The mean EW 
of the Hq, line in HD 66141 is measured to be 1 135.7 ± 18.5 mA. 
The rms of 18.5 mA corresponds to a 1.6% variations in the EW. 
Figure [6] (top panel) shows the H„ EW variations as a function 
of time and the Lomb-Scargle periodogram of the H„ EW vari- 
ations is shown in the middle panel of Figure [9] There is large 
power at the frequency of 0.001416 ± 0.000070 c d _1 with an 
FAP of less than 10~ 5 , corresponding to a period of 706.4 + 35.0 
days. Figure|7]shows the H„ EW variations phased to a period of 
706 days. 
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Fig. 7. H a EW variations for HD 66141 phased to the period of 
706 days. The same epoch of JD 2448098.0 found in the 1.25// 
(703-day flux variations) was used. 
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Fig. 8. BVS vs. RV variations for HD 66141 from December 
2003 to January 201 1. The dashed line marks a slope of 0.07. 



5.3. Line bisector variations 

The RV variations produced by the rotational modulation of the 
surface inhomogeneities should produce some changes in the 
spectral line shape, such as line asymmetry (Queloz et al. 2001). 
Thus, the variations in the shapes of spectral lines help to inter- 
pret the origin of RV variations. The difference in the bisectors 
of line widths between the top and bottom of the line profile is 
defined as the bisector velocity span (BVS). 

We measured the BVS using the least squares deconvolution 
(LSD) technique (Donati et al. 1997; Reiners & Royer 2004; 
Glazunova et al. 2008), which is calculated by the mean pro- 
file of the spectral lines. We also used the Vienna Atomic Line 
Database (VALD; Piskunov et al. 1995) to prepare the list of 
spectral lines. A total of ~ 3900 lines within the wavelength re- 
gion of 4500 - 4900 A and 6450 - 6840 A were used to con- 
struct the LSD profile, which excluded spectral regions around 
the I2 absorption region, hydrogen lines, and regions with strong 
contamination by terrestrial atmospheric lines. Then, we calcu- 
lated the BVS of the mean profile between two different central 
depth levels, 0.8 and 0.25. BVS variations as a function of time 
are exhibited in Figure [6] (bottom panel) and BVS vs. RV vari- 
ations are exhibited in Figure [8] We measured a slope of 0.07, 
which shows no correlation between BVS and the measured RV. 
The short-term accuracy of BVS measurements was measured 
by comparison of several spectra obtained in the same or dur- 
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Fig. 9. The Lomb-Scargle periodograms of the COBE/DIRBE 
satellite 1.25// flux intensity measurements, the Hipparcos pho- 
tometric measurements, the H a EW variations, and the BVS 
variations for HD 166141 (top to bottom panel), respectively. 
The vertical dashed line marks the location of the period of 480 
days and the horizontal dotted lines indicate an FAP threshold of 
1 xl0~ 3 (0.1%). Top panel - The solid line is the Lomb-Scargle 
periodogram of the COBE/DIRBE satellite 1.25// flux intensity 
measurements for 3.5 years. The arrows show the positions of 
the dominant peaks at // R = 0.00390 + 0.00008 cd _1 (Pf = 

3 IR 



0.00008 cd-' (P£ = 



256.1 ± 5.4 days) and at / 2 IR = 0.00142 : 
703.0 ± 39.4 days). The dashed line shows the periodogram of 
the residuals after removal P IR = 256 days fit to the original data. 
The dotted line shows the periodogram after removal P}£ = 703 
days fit to the original data. 



ing consecutive nights and is value of 45 m s . The long-term 
BVS scatter might be indeed larger than that of the short-term 
measurements. We found a scatter of 116 m s~' for the whole 
BVS measurements of HD 66141. In Figure [9] (bottom panel), 
the Lomb-Scargle periodogram of the BVS shows no statisti- 
cally significant peaks in the domain of interest. 
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Table 3. Orbital parameters of the best fit Keplerian orbit for HD 
66141b. 



Fig. 10. RV measurements for HD 66141 phased to the orbital 
period of 480.5 days. The solid line is the orbital solution that 
fits the data with a rms of 38.8 m s _1 . 

5.4. Orbital fit 

We find the variation fitted best with a Keplerian orbit of a period 
P = 480.5 + 0.5 days, a semi-amplitude K = 146.2 + 2.7 m s _I 
and an eccentricity e = 0.07 + 0.03. Solid line in Figure[TJshows 
the RV curve as a function of time for HD 66141 and the resid- 
uals after extracting the main frequency. The RV measurements 
phased to P — 480.5 days are shown in Figure [TUJ and they ex- 
hibit clear periodic variations at this frequency. This Keplerian 
orbit determines the minimum mass of a planetary companion m 
sin i = 6.0 + 0.3 M Jup at a distance a = 1.2 ± 0.1 AU from HD 
66141. All the orbital elements are listed in Table|3] 

The dispersion of the RV residuals is 38.8 m s _1 , which is 
significantly larger than the rms scatter of the RV standard star 
(6.7 m s -1 ) and the typical RV measurements error (~ 8 m s _1 ). 
We thus can propose that the large scatter in the residuals can be 
attributed to unresolved oscillations. Using the scaling relation 
(equation [7] from Kjeldsen & Bedding 1995) with a luminosity 
and a mass of HD66141 (Tabled, we expect the semi-amplitude 
of pulsations to be 39 m s -1 , which is almost in agreement with 
38.8 m s _1 dispersion of the residuals. 

The visual inspection of the residuals shows also a long- 
period about + 50 m s _1 variations with a time-scale comparable 
to the time span of observations. It is premature to discuss the 
origin of long-term variations until we get longer time-span ob- 
servations covering several variability cycles. For completeness, 
we would like to mention two other exoplanet-hosting candidate 
K giant stars from our survey: y l Leo (Han et al. 2010) and a Ari 
(Kim et al. 2006; Lee et al. 2011). They also exhibit long-term 
variations of the RV residuals that are superimposed with short- 
term pulsational variability well established from night-to-night 
observations. 



6. Discussion and conclusion 

From the analysis of the eight-year precise RV measurements, 
we found compelling evidence for a low-amplitude and long- 
period 480-day RV variations in a K giant HD 66141. We ex- 
amined possible origins of these RV variations. The Ca II H line 
profile, the Hipparcos photometric measurements, and the BVS 
measurements show no detectable indication of variability at this 
period. 

However, the analysis of 1 .25// COBE/DIRBE infrared flux 
measurements for HD 66141 revealed two periods of 256 or 703 



Parameter 



Value 



Period (days) 

7*periastron (ID) 

K (ms- 1 ) 
e 

0) (deg) 
Km) (M ) 
a sin i (AU) 
o-(O-C) (ms- 1 ) 



480.5 ± 0.5 
2451320.8 ± 4.8 (1999.386) 
146.2 ± 2.7 
0.07 ± 0.03 
22.1 + 3.5 
(1.545) x 10~ 7 
(6.444) x 10~ 3 
38.8 



with M* = 1.1 
m sin i (Mj up ) 
a (A\T) 



:0.1(M o ) 



6.0 : 
1.2: 



0.3 

0.1 



days that well fit the observed flux variations. Judging which 
a period is real came from the EW analysis of H a line, which 
revealed almost the same period of variations (703.0 + 39.4 
vs. 706.4 + 35.0 days). Comparison of Figure|3] and Figure|7] 
clearly shows about 0.18 phase lag between the time of maxi- 
mum H tl , EW and the maximum of the 1.25// flux. We note that 
the wavelength-to-wavelength phase lags are clearly present in 
variations of M giants and Mira variables with the optical max- 
imum appearing about 0.18 phase before the maxima at 1.25 // 
(Price et al. 2010). It was shown by Alvarez & Plez (1998) that 
optical/NIR lags in oxygen-rich Miras are likely due to strong 
titanium oxide (TiO) variability during a pulsation cycle. For the 
case of K2 III star HD 66141 the TiO lines might be formed in a 
cool spot on the surface having the temperature contrast of about 
1000 K. Thus, the temperature in the cool spot might approach to 
about 3500 K (the typical temperature of M giants). The question 
on the origin of 0.18 phase lag between H a EW and the 1.25// 
flux variations need further confirmations and investigations. 

Returning to about 703 - 706 day periods of variations in the 
NIR fluxes and optical absorptions in a spectral line we note that 
it stands far aside the RV period of 480 days. The most likely 
origin of 703 - 706 day periods is certainly the surface spot(s) 
and the rotation of HD 66141. These periods are equal to each 
other within accuracy, thus we deduce a mean value P rot = 705 
+ 53 days as the rotation period of star. With a new estimation 
of rotation period of HD 66141 and accurate radius, we can ac- 
curately re-estimated the rotational velocity as v 10t sin i = 1.5 
km s . 

In summary, the most likely cause of 480-day RV variations 
is a planet orbiting the star. HD 66141 is a planetary system with 
a 1.1 + 0.1 M giant star and a 6.0 + 0.3 M Jup planet. 
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